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Do neutrinos have nonzero masses? If they do, then these masses are very tiny, and 
can be sought only in very sensitive experiments. The most sensitive of these search 
for neutrino oscillation, a quantum interference effect which requires neutrino mass 
and leptonic mixing. In these lectures, we explain what leptonic mixing is, and 
then develop the physics of neutrino oscillation, including the general formalism 
and its application to special cases of practical interest. We also see how neutrino 
oscillation is affected by the passage of the oscillating neutrinos through matter. 



1 Introduction 

We humans, and the everyday objects around us, are made of nucleons and 
electrons, so these particles are the ones most familiar to us. However, for 
every nucleon or electron, the universe as a whole contains around a billion 
neutrinos. In addition, each person on earth is bombarded by ~ 10 14 neutrinos, 
coming from the sun, every second. Clearly, neutrinos are abundant. Thus, it 
would be very nice to know something about them. 

One of the most basic questions we can ask about neutrinos is: Do they 
have nonzero masses? Until recently, there was no hard evidence that they do. 
It was known that, at the heaviest, they are very light compared to the quarks 
and charged leptons. But for a long time, the natural theoretical prejudice 
has been that the neutrinos are not massless. One reason for this prejudice 
is that in essentially any of the grand unified theories that unify the strong, 
electromagnetic, and weak interactions, a given neutrino v belongs to a large 
multiplet F, together with at least one charged lepton £, one positively-charged 
quark q + , and one negatively-charged quark q~: 



F 



(1) 



The neutrino v is related by a symmetry operation to the other members of 
F, all of which have nonzero masses. Thus, it would be peculiar if v did not 

a To appear in the Proceedings of TASI 2000, the Theoretical Advanced Study Institute in 
elementary particle physics held in Boulder, Colorado in June 2000. 
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have a nonzero mass as well. 

If neutrinos do have nonzero masses, we must understand why they are 
nevertheless so light. Perhaps the most appealing explanation of their lightness 
is the "see-saw mechanism" Id To understand how this mechanism works, let us 
note that, unlike charged particles, neutrinos may be their own antiparticles. 
If a neutrino is identical to its antiparticle, then it consists of just two mass- 
degenerate states: one with spin up, and one with spin down. Such a neutrino 
is referred to as a Majorana neutrino. In contrast, if a neutrino is distinct from 
its antiparticle, then it plus its antiparticle form a complex consisting of four 
mass-degenerate states: the spin up and spin down neutrino, plus the spin up 
and spin down antineutrino. This collection of four states is called a Dirac 
neutrino. In the see-saw mechanism, a four-state Dirac neutrino Af D of mass 
Md gets split by "Majorana mass terms" into a pair of two-state Majorana 
neutrinos. One of these Majorana neutrinos, v u , has a small mass M v and is 
identified as one of the observed light neutrinos. The other, N M , has a large 
mass Mm characteristic of some high mass scale where new physics beyond 
the range of current particle accelerators, and responsible for neutrino mass, 
resides. Thus, N M has not been observed. The character of the breakup of 
M D into v M and N M is such that M V M^ — M|>. It is reasonable to expect 
that Mr>, the mass of the Dirac particle AT D , is of the order of Mi or q , the mass 
of a typical charged lepton £ or quark g, since the latter are Dirac particles 
too. Then M V M N = Mf m . With M im q a typical charged lepton or quark 
mass, and Mjv very big, this "see-saw relation" explains why M v is very tiny. 
Note that the see-saw mechanism predicts that each light neutrino v M is a 
two-state Majorana neutrino, identical to its antineutrino. 

2 Neutrino Oscillation 

To find out whether neutrinos really do have nonzero masses, we need an ex- 
perimental approach which can detect these masses even if they are very small. 
The most sensitive approach is the search for neutrino oscillation. Neutrino 
oscillation is a quantum interference phenomenon in which small splittings be- 
tween the masses of different neutrinos can lead to large, measurable phase 
differences between interfering quantum-mechanical amplitudes. 

To explain the physics of neutrino oscillation, we must first discuss leptonic 
"flavor" . Suppose a neutrino v is born in the VK-boson decay 

W+ £+ + v . (2) 

Here, a — e, [i or r, and is one of the positively charged leptons: £j = e + , 
= fi + , and = r + . Suppose that, without having time to change its 
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character, the neutrino v interacts in a detector immediately after its birth 
in the decay (|^), and produces a new charged lepton £7 via the reaction v 
+ target — > £J } + recoils. It is found that the "flavor" f3 of this new charged 
lepton is always the same as the flavor a of the charged lepton with which 
v was born. It follows that the neutrinos produced by the W decays (|]) to 
charged leptons of different flavors must be different objects. We take this fact 
into account by writing these decays more accurately as 

W + -> £+ + v a ; a = e,n,r . (3) 

The neutrino v ai called the neutrino of flavor a, is by definition the neutrino 
produced in leptonic W decay in association with the charged lepton of flavor 
a. As we have said, when v a interacts to create a charged lepton, the latter 
lepton is always £ a . In neutrino oscillation, a neutrino born in association with 
a charged lepton £ a of flavor a then travels for some time during which it can 
alter its character. Finally, it interacts to produce a second charged lepton £p 
with a flavor /3 different from the flavor a of the charged lepton with which 
the neutrino was born. For example, suppose a neutrino is born with a muon 
in the pion decay ir + — » Virtual W + — > fi + + v^. Suppose further that after 
traveling down a neutrino bcamlinc, this same neutrino interacts in a detector 
and produces, not another muon, but a r~ . At birth, the neutrino was a v^. 
But by the time it interacted in the detector, it had turned into a v r . One 
describes this metamorphosis by saying the neutrino oscillated from a into 
a v r . As we will see, the probability for it to change its flavor does indeed 
oscillate with the distance it travels before interacting. As we will also see, the 
oscillation in vacuum of a neutrino between different flavors requires neutrino 
mass. 

To see how neutrino mass can lead to neutrino oscillation, let us briefly 
recall the weak interactions of quarks. As we all know, there are three quarks — 
the u (up), c (charm), and t (top) quarks — which carry a positive electric charge 
Q = +2/3. In addition, there are three quarks — the d (down), s (strange), and 
b (bottom) quarks — which carry a negative electric charge Q — —1/3. Each 
of these six quarks is a particle of definite mass. As we know, the quarks are 
arranged into three families or generations, each of which contains one positive 
quark and one negative quark: 

Family : 1 2 3 



Quarks : 




However, we know experimentally that under the weak interaction, any of the 
negative quarks, d, s, or b, can absorb a positively-charged W boson and turn 
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Figure 1: Absorption of a W boson by a quark. 

into any of the positive quarks, u, c, or t. This is illustrated in Fig. ([!]). There, 
di(i — d, s, b) is one of the down-type (negative) quarks. That is, dd = d is the 
down quark, d s = s is the strange quark, and so on. Similarly, u a {a = u,c,t) 
is one of the up-type (positive) quarks, with u c = c being the charm quark, 
and so on. _ 

In the S(tandard) M(odel) of the electroweak interactionsp the VF-quark 
couplings depicted in Fig. (Q) are described by the Lagrangian density 

C udw = — % "^7V Qi cZ Ll PF+ + h.c. . (4) 

i=d,s,b 

Here, the subscript L denotes left-handed chiral projection. For instance, 
d^ s = i(l — 7s)d s is the left-handed strange-quark field. The constant g 
is the semiweak coupling constant, and V is a 3 x 3 matrix known as the quark 
mixing matrix. In the SM, V is unitary, because it is basically the matrix for 
the transformation from one basis of quantum states to another. 

The SM interaction (|]) is very well confirmed experimentally. With this es- 
tablished behavior of quarks in mind, let us now return to the leptons. Like the 
quarks, the charged leptons e, /i, and t are particles of definite mass. However, 
if leptons behave as quarks do, then the neutrinos v e ,v lll and v T of definite 
flavor are not particles of definite mass. Let us call the neutrinos which do 
have definite masses Vi, i = 1, 2, . . . , N. As far as we know, the number of 
v i: N, may exceed the number of charged leptons, three. Now, as we recall, 
the neutrino v a of definite flavor a is the neutrino state that accompanies the 
definite- mass charged lepton l a in the decay W — > l a + v a . If leptons behave 
as quarks do, this neutrino state must be a superposition of the neutrinos Vi of 
definite mass. To see this, we first note that just as any negative quark di of 
definite mass can absorb a W + and turn into any positive quark u a of definite 
mass, so it must be possible for any neutrino Vi of definite mass to absorb a 
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Figure 2: Absorption of a W boson by a neutrino. 



W~ and turn into any charged lepton l~ of definite mass. This absorption is 
illustrated in Fig. (||) . We expect that in analogy with Eq. (|]) for the VF-quark 
couplings, the SM interaction that describes the H^-lepton couplings is 

i=l,.'..,iV i 

Here, U is an N x N unitary matrix which is the leptonic analogue of the quark 
mixing matrix V. The matrix U is referred to as the leptonic mixing matrixB 
If N > 3, then only the top 3 rows of U enter in the Vl^-lepton interaction, 
Eq. (|). 

The leptonic decays of the W + are governed by the second term of C( v w, 
Eq. (g). From this term, we see that when W + — » 'V Q " , the neutrino state 
\v a ) produced in association with the specific definite-mass charged lepton ^+ 
is 

k> = £t^h>- (6) 

i 

That is, the "flavor-a" neutrino \v a ) produced together with is a coherent 
superposition of the mass-eigenstate neutrinos |z^), with coefficients which are 
elements of the leptonic mixing matrix. 

What if N is bigger than three? Suppose, for example, that N = 4. 
Then, with the elements of the bottom row of U, E/iastrow,i, w e can construct 
a neutrino state 

\ V *) = ^2 U Ltraw,i Wi) ( 7 ) 

i 

which does not couple to any of the 3 charged leptons. This state is called a 
"sterile" neutrino, which just means that it does not participate in the SM weak 
interactions. It may, however, participate in other interactions beyond the SM 
whose effects at present-day energies are too feeble to have been observed. 
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Figure 3: Creation of a neutrino with a charged lepton of flavor a followed by the interaction 
of this neutrino to produce a charged lepton of flavor /3. The "Source" is the particle whose 
decay creates the neutrino, £ a , and other, unlabelled, particles. The "Target" is the particle 
struck by the neutrino to produce tp and other, unlabelled, particles. "A" denotes an 
amplitude. 



Owing to the leptonic mixing described by Eq. (||), when the charged 
lepton of flavor a is created, the accompanying neutrino can be any of the 
z/j. Furthermore, if this Vi later interacts with some target, it can produce a 
charged lepton £p of any flavor (3. In such a sequence of events, the neutrino 
itself is an unseen intermediate state. Thus, as shown in Fig. (^), the amplitude 
for a neutrino to be born with charged lepton £ a , and then to interact and 
produce charged lepton £p, is a coherent sum over the contributions of all the 
unseen mass eigenstates Vi. 

The birth of a neutrino with charged lepton ta and its subsequent in- 
teraction to produce charged lepton A is usually described as the oscillation 
v a —> Vjj of a neutrino of flavor a into one of flavor /3 (see earlier discussion). 
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Using "A" to denote an amplitude, we see from Fig. (g) that 
A(v a — > up) — ^^[^(neutrino born with l\ is ui) x 

i 

A(vi propagates ) A(when Vi interacts it makes t^)\ . (8) 
From Civw, Eq. (g), we find that apart from irrelevant factors, 

^(neutrino born with is vi) = . (9) 

Similarly, 

A(when Vi interacts it makes £p) = Ufa . (10) 

To find the amplitude A (^..propagates), we note that in the rest frame of V{, 
where the proper time is r,p, Schrodinger's equation states that 

*^-h(ri))=M i |i/ i (T i )). (11) 
on 



Here, M» is the mass of V{. From Eq. (jll|), 

Wi(ri)) = e- lM > T ' \ Vi (0)) . (12) 

Now, for propagation over a proper time interval tj, A(i/j propagates) is just 
the amplitude for finding the original state |z^(0)) in the time-evolved \vi(ji)). 
That is, 

A(yi propagates) = (^(0)|^(r 4 )) = e~ lM ^ . (13) 

In terms of the time t and position L in the laboratory frame, the Lorentz- 
invariant phase factor exp {—iMiri) is 

e - t (E t t- Pt L) ( W ) 

Here, Ei and pi are, respectively, the energy and momentum of Vi in the 
laboratory frame. In practice, our neutrino will be highly relativistic, so if it 
was born at (t, L) — (0, 0), we will be interested in evaluating the phase factor 
( fl4|) where t w L, where it becomes 

e -<(E 4 - tt )£ _ ( 15 ) 

Suppose that the neutrino created with t a is produced with a definite 
momentum p, regardless of which v$ it happens to be. Then, if it is the 
particular mass eigenstate Vi, it has total energy 
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Ei = ^p 2 +M?np+^ , (16) 



assuming that all the masses Mi are much smaller than p. From ( |l5| ) , we then 
find that 

M 2 

A(i>i propagates) w e" 4 ^p" L . (17) 

Alternatively, suppose that our neutrino is produced with a definite energy 
E, regardless of which vi it happens to bed Then, if it is the particular mass 
eigenstate Vi, it has momentum 

Pl = ^E^-Mf^E-^l. (18) 
From @, we then find that 

M 2 

A{i>i propagates) k, e ~ l ^ L . (19) 
Since highly relativistic neutrinos have E ss p, the propagation amplitudes 



given by Eqs. ( |17| ) and (19) are approximately equal. Thus, it doesn't matter 
whether our neutrino is created with definite momentum or definite energy. 

Collecting the various factors that appear in Eq. (||) , we conclude that the 
amplitude A(y a — > vp) for a neutrino of energy E to oscillate from a v a to a 
vp while traveling a distance L is given by 

A(v a -> = £ U* ai e- iM l & . (20) 

i 

The probability P(y a — > vp) for this oscillation is then given by 
P{y a -» ^) = \A(y a -» ^) | 2 

= <5 Q(3 - 4]T SUU^UfHUajUfc) sin 2 (SM 2 A) 

+2 £ QiU^UfnU^Ufc) S in(<5M?. A) . ( 2 1) 

«>i 

Here, 5My = M 4 2 — M?, and in calculating |.A(i/ a — > ^/3)| 2 , we have used the 
unitarity constraint 

i 

Some general comments arc in order: 



1. From Eq. (21) for P(y a — > i/p), we see that if all neutrino masses vanish, 
then P(y a — > v$) — 8 a p 1 and there is no oscillation from one flavor to 
another. Neutrino flavor oscillation requires neutrino mass. 



2. The probability P(v a — ► 2/3) oscillates as a function of L/E. This is why 
the phenomenon we are discussing is called "neutrino oscillation" . 

3. From Eq. ( pfj| ) for A{v a — > 1/3), we see that the L/i? dependence of 
neutrino oscillation arises from interferences between the contributions 
of the different mass eigenstates vt. We also see that the phase of the Vi 
contribution is proportional to Mf. Thus, the interferences can give us 
information on neutrino masses. However, since these interferences can 
only reveal the relative phases of the interfering amplitudes, experiments 
on neutrino oscillation can only determine the splittings SMfj = Mf — 
M"j, and not the underlying individual neutrino masses. This fact is 
made perfectly clear by Eq. (|2l]) for P(y a — * vp). 

4. With the so-far omitted factors of Ti and c inserted, 

M^A-l^CV)^. (23) 

Thus, from Eq. ( pl| ) for the oscillation probability P(y a — ► vp), we see 
that an oscillation experiment characterized by a given value of L(km) / 
E(GeV) is sensitive to mass splittings obeying 



SM*(eV 2 ) 



L(km) 



E(GeV) 



(24) 



To be sensitive to tiny SMfj, an experiment must have large L/E. In 
Table [l], we indicate the SM 2 reach implied by Eq. ( pi[ ) for experiments 
working with neutrinos produced in various ways. 

There are basically two kinds of oscillation experiments: appearance ex- 
periments, and disappearance experiments. 

In an appearance experiment, one looks for the appearance in the neutrino 
beam of neutrinos bearing a flavor not present in the beam initially. For 
example, imagine that a beam of neutrinos is produced by the decays of 
charged pions. Such a beam consists almost entirely of muon neutrinos 
and contains no tau neutrinos. One can then look for the appearance of 
tau neutrinos, made by oscillation of the muon neutrinos, in this beam. 

In a disappearance experiment, one looks for the disappearance of some 
fraction of the neutrinos bearing a flavor which is present in the beam 
initially. For example, imagine again that a beam of neutrinos is pro- 
duced by the decays of charged pions, so that almost all the neutrinos 
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Table 1: The approximate reach in 5M 2 of experiments studying various types of neutrinos. 
Often, an experiment covers a range in L and a range in E. To construct the table, we have 
used typical values of these quantities. 



Neutrinos 
(Baseline) 


L(km) 


E{GeV) 


L(km) 
E{GeV) 


SM^eV 2 ) 
Reach 


Accelerator 


1 


1 


1 


1 


(Short Baseline) 




Reactor 
(Medium Baseline) 


1 


10^3 


10 3 


io- 3 


Accelerator 
(Long Baseline) 


10 3 


10 


10 2 


io- 2 


Atmospheric 


10 4 


1 


10 4 


io- 4 


Solar 


10 8 


10~ 3 


10 11 


10" 11 



in the beam are muon neutrinos, v^. If one knows the flux that is 
produced initially, one can look to see whether some of this initial 
flux disappears after the beam has traveled some distance, and the muon 
neutrinos have had a chance to oscillate into other flavors. 

6. Even though neutrinos can change flavor through oscillation, the total 
flux of neutrinos in a beam will be conserved so long as U is unitary. To 
see this, note that 

^PK->^) = ^|AK^^)| 2 



i j 
i,3 

= J2\U m \ 2 = l- (25) 

i 

Here, we have used Eq. ( po|) for the amplitude A(y a — > vp) and the 
unitarity relations YlpUptUpj = <5y and J2 i \U a i\ 2 = 1. The result 
^2 a P(v a — ► i/g) = 1 means that if one starts with a certain number 
n of neutrinos of flavor a, then after oscillation the number of neutrinos 
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that have oscillated away into new flavors [3 ^ a, plus the number that 
have retained the original flavor a, is still n. Note, however, that some 
of the new flavors that get populated by the oscillation might be sterile. 
If they are indeed sterile, then the number of "active" neutrinos (i.e., 
neutrinos that participate in the SM weak interactions) remaining after 
oscillation will be less than n. 



2.1 Special Cases 

Let us now apply the general formalism for neutrino oscillation in vacuum to 
several special cases of practical interest. 

The simplest special case of all is two-neutrino oscillation. This occurs 
when the SM weak interaction, Eq. (||), couples two charged leptons (say, e 
and /i) to just two neutrinos of definite mass, v\ and v 2 , and only negligibly 
to any other neutrinos of definite mass. It is then easily shown that the 2x2 
submatrix 

~U el U e2 ' 

of the mixing matrix U must be unitary all by itself. This means that the 
definite-flavor neutrinos v e and are composed exclusively of the mass eigen- 
states v\ and v 2 , and do not mix with neutrinos of any other flavor. 

From Eq. ( pfj| ) for the oscillation amplitude, we have for this two-neutrino 

case 

e +iM ?&A(v e -» = U* el U,a + U^U^ 5 ^™ . (27) 



U = 



(26) 



Since U is unitary, U^U^i + U^U^ — 0, so Eq. < \27\ ) may be rewritten as 

= - 2l e- i5M ^ & U* e2 U, 2 MSM^) . (28) 

Squaring this result, using Eq. (|2^) to take account of the requisite factors of 
h and c, we find that the probability P[v e — * v^) for a v e to oscillate into a i/„ 
is given by 

P{v e „„) = 4|[/ e2 | 2 |C/ Al2 | 2 sin 2 (1.27SM 2 (eV 2 )-|^) . (29) 

Here, we have introduced the abbreviation 8M 21 = 5M 2 . From the e *-> /x 
symmetry of the right-hand side of Eq. fl2S| ) , it is obvious that 

P(v» -» v B ) = P{v e -» . (30) 
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The probability P(y a — > v a ) that a neutrino v a of flavor a = e or /i retains 
its original flavor is given by 

P(v a -► I/ a ) = 1 - P(v a -> Upyta) 

= 1 - 4|C/ Q2 | 2 (1 - |C/ Q2 |2) s in 2 (1.27<JM 2 (eV 2 )-^^) . (31) 

To obtain this expression, we have used the conservation of probability, Eq. ( p5| ) , 
the "off-diagonal" oscillation probability, Eqs. ( p9| ) and (^p|), and the unitarity 
relation \U e2 \ 2 + |C/ m2 | 2 = 1. 

The unitarity of U, Eq. J26|), implies that it can be written in the form 



U = 



3 l ^ 1+ ^) s in6» e^+vz) COS 



(32) 



Here, 9 is an angle referred to as the leptonic mixing angle and ^1,2,3 are 
phases. From Eq. @, 4\U e2 \ 2 \U^ 2 \ 2 = sin 2 20, so that P{u e -> u^), Eq.' @, 
takes the form 

-» ^) = - v e ) = sin 2 20 sin 2 (1.27 SM 2 (eV 2 ) ) • (33) 

This is the most-commonly quoted form of the two-neutrino oscillation prob- 
ability. 

A second special case which may prove to be very relevant to the real world 
is a three-neutrino scenario in which two of the neutrino mass eigenstates are 
nearly degenerate. That is, the neutrino (Mass) 2 spectrum is as in Fig. (0), 
where 

1 5Ml I = 5M| mall « 1 5Ml I S 1 8Ml 2 \ = SM 2 ig . (34) 

All three of the charged leptons, e, fi, and r, are coupled by the SM weak 
interaction, Eq. (||), to the neutrinos ^1.2,3- 

Suppose that an oscillation experiment has L/E such that 8M^ ig L /E is 
of order unity, which implies that SM 2 ma]] L/ E <C 1. For this experiment, and 
(3 a, the oscillation amplitude of Eq. (EOh is given approximately by 

e iM i*A(u a - i/^a) S (C/*!^! + U* 2 U 02 )e i5M ^ + U* 3 U 03 . (35) 
Using the unitarity constraint of Eq. (p2), this becomes 

= _2ie i4M 3 2 a At^ 3 C^ 3 s in(<5M 3 2 2 ^) . (36) 
12 




Figure 4: A three neutrino (Mass) 2 spectrum in which the V2 — v\ splitting ^^g ma jj i s much 
smaller than the splitting <5Mj| ig between U3 and the U2 — v\ pair. The latter pair may be 
at either the bottom or the top of the spectrum. 



Taking the absolute square of this relation, using |<SAf| 2 | = ^-^ligj an d insert- 
ing the omitted factors of ft and c, we find that the v a — > vp^ a oscillation 
probability is given by 

P{v a - Vf¥a ) = 4|[/ Q3 | 2 |C/ /33 | 2 sin 2 (1.27a/l ig (eV 2 )-|^) . (37) 

To find the corresponding probability P[v a — > v a ) that a neutrino of flavor 
a retains its original flavor, we simply use the conservation of probability, 
Eq. (§): 

P(y a -^v a ) = l-^2 P(v a -► up) . (38) 

From Eq. (p^) and the unitarity relation ^2 j3 ^ a | C//33 1 2 = 1 — |C/ a 3| 2 , we then 
find that 

P{v a -,/„) = 1 - 4|l/ q3 | 2 (1 - \U a3 \ 2 ) sin 2 (1.27 5Af 2 ig (eV 2 )-i^) . (39) 

Comparing Eqs. (|37| ) and (^9|), and Eqs. ( j39"l) and (|3l|), we see that in the 
three-neutrino scenario with 5Mg ma:a L/E <C 1, the oscillation probabilities are 
the same as in the two-neutrino case, except that with three neutrinos, the 
isolated one, v 3 , plays the role played by v 2 when there are only two neutrinos. 
This strong similarity between the oscillation probabilities in the two cases is 
easy to understand. When the three-neutrino case is studied by an experiment 
for which ^Mg man L/E <C 1, the experiment cannot see the splitting between 
V\ and vi. In such an experiment, there appear to be only two neutrinos with 
distinct masses: the u\ and vi pair, which looks like a single neutrino, and 
As a result, the two neutrino oscillation probabilities holdtl 
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3 Neutrino Oscillation in Matter 

So far, we have been talking about neutrino oscillation in vacuum. However, 
some very important oscillation experiments are concerned with neutrinos that 
travel through a lot of matter before reaching the detector. These neutrinos 
include those made by nuclear reactions in the core of the sun, which traverse 
a lot of solar material on their way out of the sun towards solar neutrino 
detectors here on earth. They also include the neutrinos made in the earth's 
atmosphere by cosmic rays. These atmospheric neutrinos can be produced in 
the atmosphere on one side of the earth, and then travel through the whole 
earth before being detected in a detector on the other side. To deal with the 
solar and atmospheric neutrinos, we need to understand how passage through 
matter affects neutrino oscillation. 

To be sure, the interaction between neutrinos and matter is extremely 
feeble. Nevertheless, the coherent forward scattering of neutrinos from many 
particles in a material medium can build up a big effect on the oscillation 
amplitude. 

For both the solar and atmospheric neutrinos, it is a good approximation 
to take just two neutrinos into account. Furthermore, it is convenient to treat 
the propagation of neutrinos in matter in terms of an effective Hamiltonian. To 
set the stage for this treatment, let us first derive the Hamiltonian for travel 
through vacuum. For the sake of illustration, let us suppose that the two 
neutrino flavors that need to be considered are i/ e and v^. The most general 
time-dependent neutrino state vector \v(t)) can then be written as 

Ht)) = J2 /«WK) , (40) 

a— e,/i 

where f a (t) is the time-dependent amplitude for the neutrino to have flavor a. 
If Hv (short for Ji vacuum) is the Hamiltonian for this two-neutrino system in 
vacuum, then Schrodinger's equation for \v{t)) reads 

f ^K*)> =£ i /a(t)K> 
a 

= HvHt)) = J2.MWvW f 3) 

P 

a 

y a ) ■ (41) 



E 
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Here, (T~lv)ai3 = {va\Hy\vp)- Comparing the coefficients of \v a 
ning and end of Eq. (^Tj) , we clearly have 



d 



fe(t) 



= Tiy 



fe(t) 

M*) 



at the begin- 



(42) 



where Hy is now the 2x2 matrix with elements (Hv)ap- The Schrodinger 
equation (^2|) is completely analogous to the familiar one for a spin-1/2 particle. 
The roles of the two spin states are now being played by the two flavor states. 

Let us call the two neutrino mass eigenstates out of which v e and arc 
made v\ and i>2- To find the matrix Tiy, let us assume that our neutrino has 
a definite momentum p, so that its mass-eigenstate component vi has definite 
energy E^ given by Eq. (|l6|). That is, TivWi) = E.i\ui), and the different 
mass eigenstates like the eigenstates of any Hermitean Hamiltonian, are 
orthogonal to each other. Then, in view of Eq. (^), the elements (TCv) a /3 °f 
the vacuum Hamiltonian are given by 
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(43) 



In this expression, the two-neutrino mixing matrix U may be taken from 
Eq. (32). However, we have seen that the complex phase factors in Eq. (32) 
have no effect on the two-neutrino oscillation probabilities, Eq. (^3|). Indeed, 
it is not hard to show that when there are only two neutrinos, complex phases 
in U have no effect whatsoever on neutrino oscillation. Thus, since oscillation 
is our only concern here, we may remove the complex phase factors from the U 
of Eq. ([52j). If, in addition, we relabel the mixing angle By (short for ^vacuum), 
U becomes 

cos t)y sin ay 
— sin 9y COS By 



u 



(44) 



Inserting in Eq. (|4^) the elements U a i of this matrix and the energies Ei given 
by Eq. (|l6|), we can obtain all the (7iy) a /3- 

The matrix Ti.y can be put into a more symmetric and convenient form if 
we add to it a suitably chosen multiple (AE)I of the identity matrix /. Such 
an addition will not change the predictions of TLy for neutrino oscillation. To 
see why, we note first that the identity matrix is invariant under the unitary 
transformation that diagonalizes Tiy. Thus, adding (AE)I to Tiy in the flavor 
basis, where its elements are ( r Hv)a/3j is equivalent to adding (AE)I to Tiy 
in the mass eigenstate basis, where it is diagonal. Hence, if the eigenvalues 
of Hy are E t , i = 1, 2, those of H v + (AE)I are E, + AE, i = 1,2. That is, 
both eigenvalues are displaced by the same amount, AE. To see that such a 
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common shift of all eigenvalues does not affect neutrino oscillation, suppose a 
neutrino is born at time t — with flavor a. That is, |f(0)) = \f a )- After a 
time t, this neutrino will have evolved into the state \v(t)) given, according to 
Schrodinger's equation, by 

\ v (t)) = e- iHvt \v{Q)) 

= e-*«vt £ KM) = £ U* m e-^\v t ) . (45) 

i i 

The amplitude A(y a — ► vp) for this neutrino to have oscillated into a i/p in the 
time t is then given by 

A{v a -> ^) = (^(f)) - ]T Efte""^ . (46) 

i 

Clearly, if we add (AE)I to Tiv so that its eigenvalues are replaced by 
Ei + AE, then ^4(z/ Q — > i/g) is just multiplied by the overall phase factor 
exp[— i(A£)f]. Obviously, this phase factor has no effect on the oscillation 
probability P{v a — > !/g) = |A(i/ Q — > ^/j)| 2 - Thus, the addition of (AE)I to Tiv/ 
does not affect neutrino oscillation. 

For our purposes, the most convenient choice of AE is — [p+(M^+M|)/4p] . 
Then, from Eq. (fl3|), the new effective Hamiltonian 7iy = 7iy + (AE)I has 
the matrix elements 

(^v) a /? = g^ri^<.Ei-[p+ 1 2 ]S a p . (47) 



From Eqs. ( |44| ) and ([l6|), this gives 
w ^2 2 i 

Hv = -4E- 



— cos 28v sin 2#y 
sin 26V cos 26V 



(48) 



Here, we have used the fact that p = E, the energy of the neutrino averaged 
over its two mass-eigenstate components. We leave to the reader the instructive 
exercise of verifying that, inserted into the Schrodinger Eq. the H' v of 

Eq. j4^) does indeed lead to the usual two-neutrino oscillation probability, 

With the Hamiltonian that governs neutrino propagation through the vac- 
uum in hand, let us now ask how neutrino propagation is modified by the pres- 
ence of matter. Matter, of course, consists of electrons and nucleons. When 
passing through a sea of electrons and nucleons, a (non-sterile) neutrino can 
undergo the forward elastic scatterings depicted in Fig. <aj. Coherent for- 
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Figure 5: Forward elastic scattering of a neutrino from a particle of matter, (a) W-exchange- 
induced scattering from an electron, which is possible only for a v e . (b) Z-cxchange-induced 
scattering from an electron, proton, or neutron. This is possible for v a = v^v^, or v T . 
According to the Standard Model, the amplitude for this Z exchange is the same, for any 
given target particle, for all three active neutrino flavors. 



ward scatterings, via the pictured processes, from many particles in a material 
medium will give rise to an interaction potential energy of the neutrino in the 
medium. Since one of the reactions in Fig. (^) can occur only for electron 
neutrinos, this interaction potential energy will depend on whether the neu- 
trino is a v e or not. The interaction potential energy for a neutrino of flavor a 
must be added to the matrix element (Tl'y)aa to obtain the Hamiltonian for 
propagation of a neutrino in matterEI 

An important application of this physics is to the motion of solar neutrinos 
through solar material. The solar neutrinos are produced in the center of 
the sun by nuclear reactions such as p + p — > 2 H + e + + v e . The neutrinos 
produced by these reactions are all electron neutrinos. Let us suppose that 
the only neutrinos with which electron neutrinos mix appreciably are muon 
neutrinos, so that we have a two-neutrino system of the kind we have just been 
discussing. The solar neutrinos stream outward from the center of the sun in 
all directions, some of them eventually arriving at solar neutrino detectors here 
on earth. The passage of these solar neutrinos through solar material on their 
way out of the sun modifies their oscillation. Any neutrino which is still a 
v ei as it was at birth, can interact with solar electrons via the W exchange of 
Fig. (5a). This interaction leads to an interaction potential energy Vw{v e ) of 
an electron neutrino in the sun. This Vw(v e ) is obviously proportional to the 
Fermi coupling constant G_f, which governs the amplitude for the process in 
Fig. (5a). It is also proportional to the number of electrons per unit volume, 
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which can contribute coherently to the forward v e scattering. One can sho 
that in the Standard Model, 



a 



V w {u e ) = V2G F N e 



(49) 



This energy must be added to (7iy) ee to obtain the Hamiltonian for propaga- 
tion of a neutrino in the sun. 

In principle, the interaction energy produced by the Z exchanges of Fig. (5b) 
must also be added to H' v . However, since these Z exchanges are both flavor 
diagonal and flavor independent, their contribution to the Hamiltonian is a 
multiple of the identity matrix. As we have already seen, a contribution of this 
character does not affect neutrino oscillation. Thus, we may safely ignore it. 

In incorporating Vw{ve) m to the Hamiltonian, it is convenient, in the 
interest of symmetry, to add as well the multiple —^Vw{ve)I of the identity 
matrix. Thus, with TL' V the Hamiltonian for propagation in vacuunx-given by 
Eq. (E§) , the Hamiltonian TL Q for propagation in the sun is given bp 



G F 





-1 



AE 



- cos 29 Q 
sin 29q 



sin 26*q 
cos 28q 



In this expression, 



(50) 



D, 



sin 2 26> + (cos 29 Q - x Q ) 



and 



where 



sin 2 29 m = 



sin 2 29 v 



sin" 29 v + (cos 29 v 
2\f2G F N e E 



Xq 



5Ml 



(51) 



(52) 



(53) 



The angle 9q is the effective neutrino mixing angle in the sun when the electron 
density is N e . 

We note that H@, Eq. (|50|), has precisely the same form as H' v , Eq. (Esh. 
The only difference between these two Hamiltonians is that the parameters — 
the mixing angle and the effective neutrino (Mass) 2 splitting out in front of 
the matrix — have different values. Of course, the electron density N e is not 
a constant, but depends on the distance r from the center of the sun. Thus, 
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the parameters 9q and (5M%i)Dq are not constant either, unlike their coun- 
terparts, Ov and SM^i, in H' v . However, let us imagine for a moment that 
N e is a constant. Then Hq, like H' v , is independent of position, and must 
lead to the same oscillation probability, Eq. (|33|), as TL' V does, except for the 
substitutions 9 -> 9 Q and SM^ -> ((SMfi)!^. That is, in matter of constant 
electron density N e , Hq leads to the oscillation probability 



P(v e — > ^/j) = P(iV p — > l/ e ) 

= sin 2 26> sin 2 [1.27 5M| 1 (eV 2 )D e 



L(km) 
^(GeV)^ 



(54) 



Now let N e vary with r as it does in the real world. However, suppose 
that it varies slowly enough that the constant- N e picture we have just painted 
applies at any given radius r, but with N e (r), hence Xq (r), slowly decreasing as 
r increases. Suppose also that SM^ and E are such that x Q (r = 0) > 008 20^- 
Then, assuming that cos2#y > 0, there must be a radius r = r c somewhere 
between r = and the outer edge of the sun, where N e — > 0, such that 
Xq(t c ) = cos2(9y. From Eq. (p2|), we see that at this special radius r c , there 
is a kind of "resonance" with sin 2 29 q = 1, even if 8y is tiny. That is, mixing 
can be maximal in the sun even if it is very small in vacuo. As a result, the 
oscillation probability, which is proportional to the mixing factor sin 2 29q as 
we see in Eq. (|54|), can be very large. 

A nice picture of this enhanced probability for flavor transitions in matter 
can be gained by considering the neutrino energy eigenvalues and eigenvectors. 
If we neglect the inconsequential Z exchange contribution, and take Hq from 
Eq. (||^), then the true Hamiltonian for propagation in the sun is 



^Truo = Hq + \p ■ 



Ml 



(55) 



since the second term in this expression was subtracted from the true Hamilto- 
nian to get Hq. Now, p + (M 2 + M|)/4p = E, the energy our neutrino would 
have in vacuum, averaged over its two mass-eigenstate components. Thus, 
from Eqs. (55) and (|50|), (Egh , (|is|), and (|53|), we have in the v e — basis 



HTrue(r) 



~E 0" 


SMl 


E 


H — 

AE 



cos 29y + 2xq (r) sin 29y 
sin 29y cos 29y 



(56) 



If we continue to assume that N e (r), and consequently Xq(t), varies slowly, 
we may diagonalize this Hamiltonian for one r at a time to see how a solar 
neutrino will behave. We find from Eq. (pq) that for a given r, the energy 
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eigenvalues E± (r) are given by 



E±(r)=E 



5_M£ 
AE 



x Q (r) ± y / (x (r)-cos26»y) 2 + 



sin 2 26»t, 



(57) 



To explore the implications of these energy levels, let us suppose that SM^ 
is such that at r = 0, where N e (r) and hence x@(r) has its maximum value, 
x Q > 1. From Eq. (||), the value of G Fl the value (~ 10 26 /cc) of N e (r = 0), 
and the typical energy E (~ 1 MeV) of a solar neutrino, we find that the 
required SM^ is of order 10 -5 eV 2 . The dominant term in the energies E±(r) 
of Eq. (|57j ) will be the first one, E (~ 1 MeV). However, very interesting 
physics will result from the second term, despite the fact that this term is only 
of order 8M%JAE ~ (l(T 5 cV 2 )/lMeV - 1CT 17 MeV! 

The neutrino states which propagate in the sun without mixing signifi- 
cantly with each other are the eigenvectors of 7iTruc(Y)- To study these eigen- 
vectors, let us assume for simplicity that the vacuum mixing angle 9y is small. 
Then it quickly follows that, except in the vicinity of the special radius r c where 
i Q (r c ) = cos2#y, one of the eigenvectors of HT ruc (r), Eq. (|5rj), is essentially 
pure v ei while the other is essentially pure v^. The evolution of a neutrino 
traveling outward through the sun is then as depicted in Fig. (||). The neu- 
trino follows the trajectory indicated by the arrows. Produced by some nuclear 
process, the neutrino is born at small r as a v e . From Eq. ([56]), the eigenvector 
that is essentially v e at r — is the one with the higher energy, E + (r). Thus, 
our neutrino begins its outward journey through the sun as the eigenvector be- 
longing to the upper eigenvalue, E + . Since the eigenvectors do not cross at any 
r and do not mix appreciably, the neutrino will remain this eigenvector. How- 
ever, the flavor content of this eigenvector changes dramatically as the neutrino 
passes through the region near the radius r = r c where Xq(t) — cos2#y = 0. 
For r <C r c , the eigenvector belonging to E+(r) is essentially a v e , as one may 
see from Eq. ([36|) if one neglects the small off-diagonal terms proportional to 
sm26v, and imposes the small-r condition xq > cos20y. But for r » r c , the 
eigenvector belonging to the higher energy level E + (r) is essentially a v^, as 
one may also see from Eq. (|5^) if one neglects the off-diagonal terms and im- 
poses the large-r condition xq < cos2#y. Thus, the eigenvector corresponding 
to the higher energy eigenvalue E+(r), along which the solar neutrino travels, 
starts out as a v e at the center of the sun, but ends up as a at the outer edge 
of the sun. The solar neutrino, born a v e in the solar core, emerges from the 
rim of the sun as a f^. Furthermore, it does this with high probability even if 
the vacuum mixing angle 9y is very small, so that oscillation in vacuum would 
not have much of an effect. This very efficient conversion of solar electron 
neutrinos into neutrinos of another flavor as a result of interaction with matter 
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Figure 6: Propagation of a neutrino in the sun. The horizontal axis is linear in N e and xq. 
The solid line shows the upper eigenvalue, E+(r), for small By. The dashed lines show the 
two eigenvalues, E±(r), when there is no vacuum mixing (By = 0). A neutrino born as a iv e 
follows the path indicated by the arrows. 



is known as the Mikheyev-Smirnov- Wolfenstein (MSW) effectBO 

We turn now from solar neutrinos to atmospheric neutrinos, whose prop- 
gation through the earth entails a second important application of the physics 
of neutrinos traveling through matter. As already mentioned, an atmospheric 
neutrino can be produced in the atmosphere on one side of the earth, and then 
journey through the whole earth to be detected in a detector on the other side. 
While traveling through the earth, this neutrino will undergo interactions that 
can significantly modify its oscillation pattern. 

There is very strong evidenced that the atmospheric neutrinos born as 
muon neutrinos oscillate into neutrinos v x of another flavor, ft is known that 
v x is not a v e . It could be a i/ T , or a sterile neutrino u g , or sometimes one 
of these and sometimes the other. One way to find out what is becoming of 
the oscillating muon neutrinos is to sec whether their oscillation is afffected by 
their passage through earth-matter. The oscillation — + v s will be affected, 
but — > v T will not be. 

To see why this is so, let us first assume that the oscillation is — > v T . 
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Then the Hamiltonian He (short for 7iEarth) that describes neutrino propaga- 
tion in the earth is a 2 x 2 matrix in — v T space. Now, either a or a v r will 
interact with earth-matter via the Z exchange of Fig. (5b). This interaction 
will give rise to an interaction potential energy of the neutrino in the earth. 
However, according to the S(tandard) M(odel), the Z-exchange amplitude is 
the same for a v T as it is for a v^. Thus, in the case of — > v T oscillation, the 
contribution of neutrino-matter interaction to Kb is a multiple of the identity 
matrix. As we have alrady seen, such a contribution has no effect on oscillation. 

Now, suppose the oscillation is not — > v r but — > v s . Then He is 
a 2 x 2 matrix in — v s space. A will interact with earth-matter, as 
we have discussed, but a v s , of course, will not. Thus, the contribution of 
neutrino-matter interaction to He is now of the form 



Vz(^) 




(58) 



where Vziy^f) is the interaction potential energy of muon neutrinos produced 
by the Z exchange of Fig. (5b) . Since the matrix (^8|) is not a multiple of the 
identity, neutrino-matter interaction does affect — > v s oscillation. 

According to the SM, the forward Z-exchange amplitudes for a target e 
and a target p are equal and opposite. Thus, assuming that the earth is 
electrically neutral so that it contains an equal number of electrons and protons 
per unit volume, the e and p contributions to Vz{vp) cancel. Then Vziy^) is 
proportional to the neutron number density, N n . Taking the proportionality 
constant from the SM, we have 

Vz (^)=-^N n . (59) 

To obtain the Hamiltonian He for — > v s oscillation in the earth, we add 
to the vacuum Hamiltonian H' v of Eq. (|4jf) the contribution (|58| ) from mat- 
ter interactions, using Eq. (|5^) for Vz(v^). Of course, it must be understood 
that H'y is now to be taken as a matrix in v„ — v s space, and that the vac- 
uum (Mass) 2 splitting SM%i and mixing angle 9y in H'y are now different 
parameters than they were when we obtained from H'y the Hamiltonian Hq 
for neutrino propagation in the sun. The Quantities SM^i and 6y are now new 
parameters appropriate to the vacuum oscillation of atmospheric, rather than 
solar, neutrinos. To obtain a more symmetrical and convenient He from H'y, 
we also add — hVz(Vfj,)I, a multiple of the identity which will not affect the 
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implications of He for oscillation. The result is 



He 



(SM^De 



Here, 



-cos 26 e sin20# 
sin 26 e cos 29 e 



(60) 



D E = \fsm 2 20 v + (cos28 v -x E ) 2 , (61) 

and 

. sin 2 29v ,„„s 

sin 2 28 E = — o , 62 

sin 2 20 v + (cos2<V - x E ) 2 

where 



V2G F N n E 



x e = 7TJT • (63) 

As a rough approximation, we may take the neutron density N n to be 
constant throughout the earth. Then, xe is also a constant, and the Hamilto- 
nian He of Eq. (|6(]) is identical to the vacuum Hamiltonian H' v of Eq. (f|8j), 
except that the constant 8M 2 X is replaced by the constant (8M 2 1 )De 1 and the 
constant dy by the constant 0e- Thus, from the fact that H' v leads to the 
vacuum oscillation probability of Eq. (|33] ) (with v e — > replaced by — » i/ s 
for the present application), we immediately conclude that the He of Eq. (60) 
leads to the oscillation probability 

P(p„ - v s ) = sin 2 28 E sm 2 [l.275MUeW 2 )D E ^0-] . (64) 

As we see from Eq. (p^), xe, which is a measure of the influence of matter 
effects on atmospheric neutrinos, grows with energy E. In a mimient we will see 
that for E ~ 1 GcV, matter effects are negligible. Fits to data! 2 ] on atmospheric 
neutrinos with roughly this energy have led to the conclusion that atmospheric 
neutrino oscillation involves a neutrino (Mass) 2 splitting 8M\ tmoa given by 

6M 2 Atmos ~ 3 x 10- 3 eV 2 , (65) 

and a neutrino mixing angle #Atmos given by 

sin 2 26> A tmos ~ 1 • (66) 

That is, the mixing when matter effects are negligible is very large, and perhaps 
maximal. The quantities SM 2 trnr , s and #Atmos are to be taken, respectively, for 
8M%i and 9y in Eqs. ( |60| ) - (Q) to find the implications of those equations for 
v., — > within the earth. 
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Since atmospheric neutrino oscillation involves maximal mixing when mat- 
ter effects are negligible, the matter effects cannot possibly enhance the oscil- 
lation, but can only suppress it. From Eq. ([32]), we see that if, as observed, 
sin 2 29y — 1, then matter effects will lead to a smaller effective mixing sin 2 29 e 
in the earth, given by 

sin 2 20£ = — Up . (67) 
l + x E 

As is clear in Eq. d&jj), this will result in a smaller oscillation probability 
P{ v ii ~ * v s ) than one would have in vacuum, where sin 2 29e is replaced by 
sin 2 26»y (-1). 

Since xe grows with energy, the degree to which matter effects suppress 
— ► v s grows as well. From Eq. (p3) for xe, the known values of Gf and N n , 
and the value ( |65[ ) required for 5M£ 1 by the data, we find that when 
E ~ 1 GeV. Thus, matter effects are indeed negligible at this energy, so it is 
legitimate to determinaij the vacuum parameters SMj^^^ and sin 2 2#Atmos by 
analyzing the ~ 1 GeV data neglecting matter effects. However, at sufficiently 
large E, the matter-induced suppression of — > v s will obviously be signif- 
icant. From Eq. (|3|), we find that sin 2 29 E is below 1/2 when E > 20 GeV. 
The consequenLsuppression of oscillation at these energies has been looked for, 
and is not scenEj This absence of suppression is a powerful part of the evidence 
that the neutrinos into which the atmospheric muon neutrinos oscillate are not 
sterile neutrinos, or at least not solely sterile neutrinosliij 

4 Conclusion 

Evidence has been reported that the solar neutrinos, the atmospheric neu- 
trinos, the accelerator-generated neutrinos studied by the Liquid Scintillator 
Neutrino Detector (LSND) experiment at Los Alamos, and the accelerator- 
generated neutrinos studied by the K2K experiment in Japan, actually do 
oscillate. Some of this evidence is very strong. The neutrino oscillation expefc 
iments, present and future, are discussed in this Volume by John WilkersonliJ 
In these lectures, we have tried to explain the basic physics that underlies 
neutrino oscillation, and that is invoked to understand the oscillation experi- 
ments. As we have seen, neutrino oscillation implies neutrino mass and mixing. 
Thus, given the compelling evidence that at least some neutrinos do oscillate, 
we now know that neutrinos almost certainly have nonzero masses and mix. 
This knowledge rasises a number of questions about the neutrinos: 

• How many neutrino flavors, including both interacting and possible ster- 
ile flavors, are there? Equivalently, how many neutrino mass eigenstates 
are there? 
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• What are the masses, Mj, of the mass eigenstates v{t 

• Is the antiparticle VI of a given mass eigenstate the same particle as 

or a different particle? 

• What are the sizes and phases of the elements U a i of the leptonic mixing 
matrix? Equivalently, what are the mixing angles and complex phase 
factors in terms of which U may be described? Do complex phase factors 
in U lead to CP violation in neutrino behavior? 

• What are the electromagnetic properties of neutrinos? In particular, 
what are their dipole moments? 

• What are the lifetimes of the neutrinos? Into what do they decay? 

• What is the physics that gives rise to the masses, the mixings, and the 
other properties of the neutrinos? 

Seeking the answers to these and other questions about the neutrinos will 
be an exciting adventure for years to come. 
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